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We report an observation of new bottom baryons produced in pp collisions at the Tevatron. Using 
1.1 fb _1 of data collected by the CDF II detector, we observe four A^* resonances in the fully 
reconstructed decay mode A° — > A+7r~, where A+ — > pK~iv + . We interpret these states as the 
S^*' )± baryons and measure the following masses: 

m E+ = 5807.812.2 tat -) ± 1-7 (syst.) MeV/c 2 
m s _ = 5815.2 ± 1.0 (stat.) ± 1.7 (syst.) MeV/c 2 
m(EJ) - m(S 6 ) = 21.2±l° (stat.)+°'* (syst.) MeV/c 2 

PACS numbers: 14.20.Mr, 13.30.Eg 
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Recently the CDF II detector at the Fermilab Teva- 
tron has accumulated the world's largest sample of fully 
reconstructed A® baryons, which consist of the u, d, and 
b quarks, with 3180 ± 60 (stat.) Ag -> A+tt~ candi- 
dates. This is made possible by the large bb production 
cross-section in pp collisions at y/s = 1.96 TeV and the 
ability of the CDF II experiment to select fully hadronic 
decays of b hadrons with a specialized trigger system. 
In this Letter, we present an observation of four A^tt^ 
resonances, where A° — > A+7r~ and A+ — > pK~ir + , us- 
ing 1.1 fb _1 of data. The A®ir states are interpreted as 
the lowest-lying charged E& baryons and will be labeled 
££ ' . The symbol £ h refers to £^ , while ££ refers to 
E^. Any reference to a specific charge state implies the 
antiparticle state as well. 
(*)+ 

The Ej baryons contain one b and two u quarks, 
while the E^ baryons contain one b and two d quarks; 
these states are expected to exist but have not been ob- 
served. Baryons containing one bottom quark and two 
light quarks can be described by heavy quark effective 
theory (HQET) fl. In HQET a bottom baryon consists 
of a 6 quark acting as a static source of the color field 
surrounded by a diquark system comprised of the two 
light quarks. In the lowest-lying E^ states, the diquark 
system has strong isospin 1=1 and J p = 1 + , which 
couple to the heavy quark spin and result in a doublet 
of baryons with J p = ± + (E 6 ) and J p = | + (£*). This 
doublet is degenerate for infinite b quark mass. As the b 
quark mass is finite, there is a hyperfine mass splitting 
between the | and k states. There is also an isospin 

( -|- \ I 

mass splitting between the E^ and E^ states. 

Predictions for the E^ masses come from non- 
relativistic and relativistic potential quark models Q, 
1/N C expansion [3j, quark models in the HQET approx- 
imation 4], sum rules l5J, and lattice quantum chromo- 
dynamics calculations @ . On the basis of @, H, 0, IE 0] , 
we expect m(£&) - m{A° b ) - 180 - 210 MeV/c 2 , m(E*) - 
m(E 6 ) - 10 - 40 MeV/c 2 , and m(E^ 3_ ) - m(Ej,* )+ ) ~ 
5 — 7 MeV/c 2 . The difference between the isospin mass 
splittings of the ££ and E^ multiplets is predicted to 
be [m(E*+) - m(E*-)] - [m(E+) - m(Ej )] = 0.40 ± 
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0.07 MeV/c 2 • The natural width of E^ baryons is ex- 
pected to be dominated by the P-wave one pion transition 
E^ — > A^ir, whose partial width depends on the avail- 
able phase space and the pion coupling to a constituent 
quark. Using an HQET prediction [8|] , the natural widths 
for the expected E> ' masses are r(£±) « 7 MeV/c 2 and 
T(E* ± ) w 13 MeV/c 2 . 

The CDF II detector is described in detail else- 
where Q . Its components and capabilities most relevant 
to this analysis are the tracking system [l(| and a dis- 
placed track trigger which is employed to select bottom 
and charmed hadrons [ill ]. 

In reconstructing the decays A b — > A+7r~ and A+ — > 
pK~ir + , the proton from the A+ decay and the 7r~ from 
the A° decay both must have pt > 2 GeV/c [l3 |. while 
the K~ and ir + candidates have pr > 0.5 GeV/c. We 
also require Pt(p) > _Pt(tt + ) to suppress A+ combinato- 
rial background. No particle identification is used in this 
analysis. All particle hypotheses consistent with the can- 
didate decay structure are considered. In a 3-D kinematic 
fit, the A+ daughter tracks are constrained to originate 
from a single point. The A+ candidate is constrained to 
the known A+ mass, and the A+ momentum vector is ex- 
trapolated to intersect the tt~ momentum vector to form 
the A b vertex. The probability of the 3-D A b kinematic 
vertex fit must exceed 0.1%, and the A+ and A b must 
have pt greater than 4.5 and 6.0 GeV/c, respectively. 
To suppress prompt backgrounds from the primary inter- 
action, we make the following decay time requirements: 
ct(AjJ) > 250 fim and its significance ct(A b )/a ct > 10. 
We define ct(A° b ) = L xy (A°)TO A oc/p T (A°) as the A° 
proper time, where L xy (A b ) is defined as the projec- 
tion, onto pr(Aj}), of the vector connecting the primary 
vertex to the A® vertex in the transverse plane. We 
use a primary vertex determined event-by-event when 
computing this vertex displacement. To reduce com- 
binatorial backgrounds and partially reconstructed de- 
cays, we also require |d (A°)| < 80 fim, where d (A b ) 
is the impact parameter of the momentum vector of 
the A® candidate with respect to the primary vertex. 
To suppress the contributions from B° — > D + tt~ de- 
cays, where D + — > K~ir + T{ + , we require m(pK~Tr + ) 
to be within 16 MeV/c 2 of the known A+ mass 13], 
and ci(A+) G [-70,200] fim. We define ct(A+) = 
L xy (A+)m A +c/pT(A+) as the A+ proper time, where 
-L xy (A+) is defined analogously to L xy (A b ) but computed 
with respect to the A° vertex. 

The invariant mass distribution of A+7r~ candidates is 
shown in Figure [T] overlaid with a binned maximum like- 
lihood fit. A clear A° — > A+7T - signal is observed at the 
expected A b mass. The invariant mass distribution is de- 
scribed by several components: the A b — * A+7r~ signal, 
a combinatorial background, partially and fully recon- 
structed B mesons which pass the A+7r~ selection cri- 
teria, partially reconstructed A® decays, and fully recon- 
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A°7T, 



A+7T , and A+ — » 
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arc 



allowed. We optimize the T, b selection criteria by 



Total Fit 

Partially-reconstructed A b 

Fuiiy-reconstmcted b maximizing c(Smc ) /VB, where e(S'Mc) i s the efficiency 

Partially-reconstructed B v-i(*) • • rr-i • i t-> 

combinatorial 01 tne ^6 signal measured m the ML simulation and B 

Data is the number of background events in the signal region 

estimated from the upper and lower sideband regions of 
Q G [0,30] MeV/c 2 and Q G [100,500] MeV/c 2 . These 
sideband regions are parameterized by a power law 
multiplied by an exponential. We combine the A®tt~ and 
A b n + subsamples to optimize cuts on the px of the 
candidate, the impact parameter significance |do/0Vj o | °f 
FIG. 1: Fit to the invariant mass of A° -» A+7r" candidates. the 7r Sb track, and the cos#* of the 7T Si; track, where 
Curves for fully reconstructed A° decays such as A° — > A+7r~ 0* is defined as the angle between the momentum of 
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and A° — > A^ K are not indicated on the figure. The A° £h e ^ m ^ e £j 



(*) 



signal region, m(A+7r~) G [5.565,5.670] GeV/c 2 , consists 
primarily of A° baryons, with some contamination from B 
mesons and combinatorial events. The discrepancies between 
the fit and data below the AjJ signal region are due to in- 
complete knowledge of the branching ratios of the decays in 



this region and are included in the T, b 
systematics. 



M 



background model 



structed A° decays other than A+7r~ (e.g. A® — > A+AT~). 
The combinatorial background is modeled with an expo- 
nentially decreasing function. All other components are 
represented in the fit by fixed shapes derived from Monte 
Carlo (MC) simulations 14 , 1]| . The normalizations are 



constrained by Gaussian terms to branching ratios that 
are either measured (for B meson decays) or theoretical 
predictions (for A b decays). In the fit, the A b components 
are normalized relative to the A° — > A+7r~ signal. To 
normalize the B meson components, we explicitly recon- 
struct a B° — ► (K~ir + Tr + )Tr~ signal in the A+7r~ sample 
by replacing the proton mass hypothesis with the pion 
mass hypothesis. The fit to the invariant A+7r _ mass 
distribution results in 3180 ± 60 (stat.) Ag -» A+tt" 
candidates. 



The reconstruction of £^ proceeds by combining 
A° candidates in the A b signal region, m(A+7r~) G 
[5.565,5.670] GeV/c 2 , with all remaining high quality 
tracks. A pion mass hypothesis is used when comput- 
ing the invariant mass of the EJ, candidate. We search 
for narrow resonances in the mass difference distribution 
of Q = to(A°7t) — m(A b ) — m 7ll where m(A b y ) is the recon- 
structed A+7T - mass. The Ej; candidates are divided 
into two subsamples using the charge of the pion from 



E 



(*) 



b decay, denoted by 7T£ b : in the A b ir~ subsample the 
w-£ b has the same charge as the pion from A b while in the 
A°7T + subsample the has the opposite charge as the 
pion from A®. 



The 
[30,100 
in 0, 

PYTHIA 



.(*) 



signal 

,2 



)1 MeV/c 

is, a m 



region, defined as Q G 
is motivated by the predictions 
The signal is modeled by the 
event generator where only the decays 



total momentum in the lab frame 



rest frame and the direction of the 
The maximum 

of e(S M c)/VB is realized for pt(E 6 ) > 9.5 GeV/c, 
\d /a do \ < 3.0, and costf* > -0.35. 

(*) 

In the E^ search, the dominant background is from 
the combination of prompt A? baryons with extra tracks 
produced in the hadronization of the b quark. The 
remaining backgrounds arc from the combination of 
hadronization tracks with B mesons reconstructed as 
A® baryons, and from combinatorial background events. 
The percentage of each background component in the 
A® signal region, computed from the A^ mass fit, is 
(89.5 ± 1.7)% A° baryons, (7.2 ± 0.6)% B mesons, 
and (3.3 ± 0.1)% combinatorial events. Other back- 
grounds such as 5-track decays of B + mesons are neg- 
ligible, as confirmed in inclusive single b hadron sim- 
ulations 14, [TBj]. The high mass region above the 
A® — ► A+7T~ signal in Figure [TJ determines the combina- 
torial background shape. Reconstructing B° — > D + ir~ 
data as A b — ► A+7r~ gives the B hadronization back- 
ground shape. The A® hadronization background shape 
is obtained from a A b — > A+7r~ pythia simulation. 
The events in this simulation are reweighted so that 
the pr(A°) distribution agrees with data. As the sim- 
ulation has fewer low momentum tracks around the A® 
than found in data, the simulated events are further 
reweighted until the pt spectrum of tracks around the 
A° is consistent with data. The background shapes are 
parameterized by a power law multiplied by an exponen- 
tial, and the normalizations are fixed from the percentage 
of that background component in the A® signal region. 
The total background shown in Figure [2] (inset) is com- 
patible with the Q sidebands and the background shape 
and normalization are fixed components of the T, b fit. 

In the Q signal region we observe an excess of events 
over the total background as shown in Figure [2] The 
excess in the A b ir~ subsample is 118 over 288 expected 
background candidates. In the A b ir + subsample the ex- 
cess is 91 over 313 expected background candidates. 

We perform a simultaneous unbinned maximum likeli- 
hood fit to the A b ir~ and A°7r + subsamples for a signal 
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50 100 150 200 

Q = m(A°ji) - m(A°) - m n (MeV/c 2 ) 

FIG. 2: The E^*' fit to the A£tt + and A£tt~ subsamples. The 
top plot shows the A°n + subsample, which contains E^ + , 
while the bottom plot shows the A®tv~ subsample, which con- 
The insets show the expected background plotted 



(*)- 



tains E u 

on the data for Q £ [0, 500] MeV/c 2 , while the signal fit is 
shown on a reduced range of Q G [0, 200] MeV/c 2 . 



from each expected Tr b state plus the background, re- 
ferred to as the "four signal hypothesis." Each signal con- 
sists of a non-relativistic Breit-Wigner distribution con- 
voluted with two Gaussian distributions describing the 
detector resolution, with a dominant narrow core of an 
1.2 MeV/c 2 width and a small broad component of a 
3 MeV/c 2 width for the tails. The natural width of each 
Breit-Wigner distribution is computed from the central 
Q value [8J. The expected difference of the isospin mass 
splittings within the ££ and £;, multiplets is below our 
sensitivity with this sample of data. Consequently, we 
constrain m(Y>* b + ) — to(E^) = m(T,* b ~) — m(E^) = As* . 
The four Ej, signal fit to data, which has a fit probability 
of 76% in the range Q € [0,200] MeV/c 2 , is shown in 
Figure [U 

Systematic uncertainties on the mass difference and 
yield measurements fall into three categories: mass scale, 



£i background model, and X^* -1 signal parameteriza- 
tion. The systematic uncertainty on the mass scale is 
determined by the discrepancies of the CDF II mea- 
sured Q values of the D*. E e , and A* hadrons from the 
world average Q values 131 ] . The Q value dependence 



(*) 



of this systematic uncertainty is modeled with a linear 
function, which is used to extrapolate the mass scale 
uncertainty for each Tr b Q value. This is the largest 
systematic uncertainty for the mass difference measure- 



ments, ranging from 0.1 to 0.3 MeV/c 2 . The systematic 
effects related to assumptions made on the E^ back- 
ground model are: the sample composition of the A b sig- 
nal region, the normalization and functional form of the 
A b hadronization background taken from a pythia sim- 
ulation, and our limited knowledge of the shape of the A b 
hadronization background (the largest systematic uncer- 
tainty on the yield measurements, ranging from 2 to 15 
events). The systematic effects related to assumptions 
made on the E^ signal parameterization are: underesti- 
mation of the detector resolution, the uncertainty in the 
natural width prediction from [8[ , and the constraint that 
m(E h *+) - m(£+) = m(S*-) " m(Ej ). 

The significance of the signal is evaluated using the 
likelihood ratio, LR = L/L a \ t , where L is the likelihood 
of the four signal hypothesis and L a \ t is the likelihood 
of an alternative hypothesis [13] • We study the alter- 
nate hypotheses of no signal, two E/, states (one per A b ir 

(*) 

charge combination), and three E^ states, performed 
by eliminating one of the states in the four signal hy- 
pothesis. Systematic variations are included in the fit 
as nuisance parameters over which the likelihood is inte- 
grated. The resulting likelihood ratios are given in Table 
HJ To assess the significance of the signal, we repeat the 
four signal hypothesis fit on samples randomly generated 
from alternate signal hypotheses. In 12 million back- 
ground samples, none had a LR equivalent or greater 
than the one found in data. We evaluate the probabil- 
ity for background only to produce four signals of this 
or greater significance to be less than 8.3 x 10~ 8 , cor- 
responding to a significance of greater than 5.2 o. The 
probabilities for each of the alternate hypotheses to pro- 
duce the observed signal structure is also given in Ta- 
ble H] The final results for the E^ measurement are 
quoted in Table |TTJ Using the CDF II measurement of 
m A o = 5619.7 ± 1.2 (stat.) ± 1.2 (syst.) MeV/c 2 [3, we 
find the absolute masses of the E& states given in Table 
HTl The systematic uncertainties on the absolute E& mass 
values are dominated by the total A b mass uncertainty. 



TABLE I: Likelihood ratios (LR) in favor of the four signal 
hypothesis over alternative hypotheses. Also shown is the 
probability for each hypothesis to produce the observed data 
(p- value), calculated using the LR as a test statistic on ran- 
domly generated samples. The final column gives the equiv- 
alent standard deviations from the normal distribution. 



Hypothesis 


LR 


p- value 


Significance (a) 


No Signal 


2.6 x 10 18 


< 8.3 x 10~ s 


> 5.2 


Two Eb States 


4.4 x 10 6 


9.2 x 1Q- 5 


3.7 


No Ej Signal 


1.2 x 10 5 


3.2 x 10~ 4 


3.4 


No E+ Signal 


49 


9.0 x 10~ 3 


2.4 


No E*~ Signal 4.9 x 10 4 


6.4 x 10~ 4 


3.2 


No E*+ Signal 


8.1 x 10 4 


6.0 x 10~ 4 


3.2 



In summary, using a sample of 3180 ± 60 (stat.) A" 
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TABLE II: Final results for the E& measurement. The first 
uncertainty is statistical and the second is systematic. The 
absolute Ef, mass values are calculated using a CDF II mea- 
surement of the mass which contributes to the sys- 
tematic uncertainty. 



State 


Yield 


Q or 


A s; (MeV/c^) 


Mass (MeV/c^) 


y + 


oo+13+b 

JZ -12-3 


Q s + = 


,0 r+2. () + (). 2 

*0.0_ 2 2-0.3 


5807.8;2;2 ±L7 


^b 


OJ -14-4 


Q s - = 


: 55.9 ± 1.0 ± 0.2 5815.2 ± 1.0 ± 1.7 


^b 

n~ 


77 +17+10 
' '-16-6 

fi q+i8+i6 

UJ -17-5 




91 9+2.0+0.4 
zi - z -1.9-0.3 


5829.0±i;l±ii 
5836.4 ±2.01^ 



Aj 



candidates reconstructed in 1.1 fb of CDF II 



data, we search for resonant A°7r ± states. We observe a 
signal of four states whose masses and widths are con- 
sistent with those expected for the lowest-lying charged 
baryons. This result represents the first observation 

of the E^ baryons. 
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